Dark matter subhalos and the dwarf satellites of the Milky Way

The Via Lactea simulation of the dark matter halo of the Milky Way predicts the existence of many thousands of bound subhalos distributed approximately with equal mass per decade of mass. Here we show that: a) a similar steeply rising subhalo mass function is also present at redshift 0.5 in an elliptical-sized halo simulated with comparable resolution in a different cosmology. Compared to Via Lactea, this run produces nearly a factor of two more subhalos with large circular velocities; b) the fraction of Via Lactea mass brought in by subhalos that have a surviving bound remnant today with present-day peak circular velocity Vmax>2 km/s (>10 km/s) is 45% (30%); c) because of tidal mass loss, the number of subhalos surviving today that reached a peak circular velocity of >10 km/s throughout their lifetime exceeds half a thousand, five times larger than their present-day abundance and more than twenty times larger than the number of known satellites of the Milky Way; e) unless the circular velocity profiles of Galactic satellites peak at values significantly higher that expected from the stellar line-of-sight velocity dispersion, only about one in five subhalos with Vmax>20 km/s today must be housing a luminous dwarf; f) small dark matter clumps appear to be relatively inefficient at forming stars even well beyond the virial radius; g) the observed Milky Way satellites appear to follow the overall dark matter distribution of Via Lactea, while the largest simulated subhalos today are found preferentially at larger radii; h) subhalos have central densities that increase with Vmax and reach 0.1-0.3 Msun/pc3 comparable to the central densities inferred in dwarf spheroidals with core radii >250 pc.Comment: 14 pages, 8 figures, ApJ in press. A few typos correcte

On the age-radius relation and orbital history of cluster galaxies

We explore the region of influence of a galaxy cluster using numerical simulations of cold dark matter halos. Many of the observed galaxies in a cluster are expected to be infalling for the first time. Half of the halos at distances of one to two virial radii today have previously orbited through the cluster, most of them have even passed through the dense inner regions of the cluster. Some halos at distances of up to three times the virial radius have also passed through the cluster core. We do not find a significant correlation of ``infall age'' versus present day position for substructures and the scatter at a given position is very large. This relation may be much more significant if we could resolve the physically overmerged galaxies in the central region.Comment: To appear in the proceedings of IAU Colloquium 195: "Outskirts of galaxy clusters: intense life in the suburbs", Torino, Italy, March 12-16, 200

Black Holes in our Galactic Halo: Compatibility with FGST and PAMELA Data and Constraints on the First Stars

10 to 10^5 solar mass black holes with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today are examined in light of recent data from the Fermi Gamma-Ray Space Telescope (FGST). The dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation. We examine the signatures of annihilations into gamma-rays, electrons and positrons, and neutrinos. We find that some significant fraction of the point sources detected by FGST might be due to dark matter annihilation near black holes in our Galaxy. We obtain limits on the properties of dark matter annihilations in the spikes using the information in the FGST First Source Catalog as well as the diffuse gamma-ray flux measured by FGST. We determine the maximum fraction of high redshift minihalos that could have hosted the formation of the first generation of stars and, subsequently, their black hole remnants. The strength of the limits depends on the choice of annihilation channel and black hole mass; limits are strongest for the heaviest black holes and annhilation to $b \bar{b}$ and $W^+W^-$ final states. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted and thermonuclear burning runs its course; thus FGST observations may be used to constrain the properties of Dark Stars. Additionally, we comment on the excess positron flux found by PAMELA and its possible interpretation in terms of dark matter annihilation around these black hole spikes.Comment: 34 pages, 11 figures. v2: typos corrected, references added. v3: updated to match published versio

Radial distribution and strong lensing statistics of satellite galaxies and substructure using high resolution LCDM hydrodynamical simulations

We analyse the number density and radial distribution of substructures and satellite galaxies using cosmological simulations that follow the gas dynamics of a baryonic component, including shock heating, radiative cooling and star formation within the hierarchical concordance LCDM model. We find that the dissipation of the baryons greatly enhances the survival of subhaloes, expecially in the galaxy core, resulting in a radial distribution of satellite galaxies that closely follows the overall mass distribution in the inner part of the halo. Hydrodynamical simulations are necessary to resolve the adiabatic contraction and dense cores of galaxies, resulting in a total number of satellites a factor of two larger than found in pure dark matter simulation, in good agreement with the observed spatial distribution of satellite galaxies within galaxies and clusters. Convergence tests show that the cored distribution found by previous authors in pure N-body simulations was due to physical overmerging of dark matter only structures. We proceed to use a ray-shooting technique in order to study the impact of these additional substructures on the number of violations of the cusp caustic magnification relation. We develop a new approach to try to disentangle the effect of substructures from the intrinsic discreteness of N-Body simulations. Even with the increased number of substructures in the centres of galaxies, we are not able to reproduce the observed high numbers of discrepancies observed in the flux ratios of multiply lensed quasars.Comment: 11 pages, 15 figures, comparison with previous works updated, one more plot added, minor changes to match the accepted version by MNRA

Dark matter substructure and gamma-ray annihilation in the Milky Way halo

We present initial results from ``Via Lactea'', the highest resolution simulation to date of Galactic CDM substructure. It follows the formation of a Milky Way-size halo with Mvir=1.8x10^12 Msun in a WMAP 3-year cosmology, using 234 million particles. Over 10,000 subhalos can be identified at z=0: Their cumulative mass function is well-fit by N(>Msub)= 0.0064 (Msub/Mvir)^(-1) down to Msun=4x10^6 Msun. The total mass fraction in subhalos is 5.3%, while the fraction of surface mass density in substructure within a projected distance of 10 kpc from the halo center is 0.3%. Because of the significant contribution from the smallest resolved subhalos, these fractions have not converged yet. Sub-substructure is apparent in all the larger satellites, and a few dark matter lumps are resolved even in the solar vicinity. The number of dark satellites with peak circular velocities above 10 km/s (5 km/s) is 124 (812): of these, 5 (26) are found within 0.1 Rvir, a region that appeared practically smooth in previous simulations. The neutralino self-annihilation gamma-ray emission from dark matter clumps is approximately constant per subhalo mass decade. Therefore, while in our run the contribution of substructure to the gamma-ray luminosity of the Galactic halo amounts to only 40% of the total spherically-averaged smooth signal, we expect this fraction to grow significantly as resolution is increased further. An all-sky map of the expected annihilation gamma-ray flux reaching a fiducial observer at 8 kpc from the Galactic center shows that at the current resolution a small number of subhalos start to be bright enough to be visible against the background from the smooth density field surrounding the observer.Comment: 15 pages, 10 figures. ApJ accepted. Movies, images and a version with higher resolution figures are available at http://www.ucolick.org/~diemand/vl

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.Comment: 15 pages, 5 figures, accepted for publication by the Astrophysical Journa

Two-body relaxation in cold dark matter simulations

N-body simulations of the hierarchical formation of cosmic structures suffer from the problem that the first objects to form always contain just a few particles. Although relaxation is not an issue for virialized objects containing millions of particles, collisional processes will always dominate within the first structures that collapse. First we quantify how the relaxation varies with resolution, softening and radius within isolated equilibrium and non-equilibrium cuspy haloes. We then attempt to determine how this numerical effect propagates through a merging hierarchy by measuring the local relaxation rates of each particle throughout the hierarchical formation of a dark matter halo. The central few per cent of the final structures - a region that one might naively think is well resolved at the final time since the haloes contain ≈106 particles - suffer from high degrees of relaxation. It is not clear how to interpret the effects of the accumulated relaxation rate, but we argue that it describes a region within which one should be careful about trusting the numerical results. Substructure haloes are most affected by relaxation since they contain few particles at a constant energy for the entire simulation. We show that relaxation will flatten a cusp in just a few mean relaxation times of a halo. We explore the effect of resolution on the degree of relaxation, and we find that increasing N slowly reduces the degree of relaxation ∝N−0.25 rather than proportional to N as expected from the collisionless Boltzmann equation. Simulated with the same relative mass resolution (i.e. equal numbers of particles), cluster mass objects suffer significantly more relaxation than galaxy mass objects since they form relatively late and therefore more of the particles spend more time in small-N haloe